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Abstract 

The effects of retardation in the Coulomb excitation of radioactive nuclei in intermediate energy 
collisions (Ei a f, ~ 100 MeV/nucleon) are investigated. We show that the excitation cross sections 
of low-lying states in Be, > °' S and 44 > 46 Ar projectiles incident on gold and lead targets are 
modified by as much as 20% due to these effects. The angular distributions of decaying gamma- 
rays are also appreciably modified. 
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I. INTRODUCTION 



The excitation of a nucleus by means of the electromagnetic interaction with another 
nucleus is known as Coulomb excitation. Since the interaction is proportional to the charge 
Z of the nucleus, Coulomb excitation is especially useful in the collision of heavy ions, with 
cross sections proportional to Z 2 . Pure Coulomb excitation is assured if the bombarding 
energy is sufficiently below the Coulomb barrier. In this case the ions follow a Rutherford 
trajectory and never come so close so that their nuclear matter overlaps. This mechanism 
has been used for many years to study the electromagnetic properties of low-lying nuclear 
states 

The probability for Coulomb excitation of a nuclear state | /) from an initial state | i) is 
large if the transition time t/j = h/(Ef — Ei) = 1/ujfi is greater than the interaction time 
t C oii — a o/v, in a heavy ion collision with closest approach distance and projectile velocity 
v. That is, the cross section for Coulomb excitation is large if the adiabacity parameter 
satisfies the condition 

This adiabatic cut-off limits the possible excitation energies below 1-2 MeV in sub-barrier 
collisions. 

A possible way to overcome this limitation, and to excite high-lying states, would be 
the use of higher projectile energies. In this case, the closest approach distance, at which 
the nuclei still interact only electromagnetically, is of order of the sum of the nuclear radii, 
R = Rp + Rt, where P refers to the projectile and T to the target. For very high energies 
one has also to take into account the Lorentz contraction of the interaction time by means 
of the Lorentz factor 7 = (1 — t> 2 /c 2 ) -1 / 2 , with c being the speed of light. For such collisions 
the adiabacity condition, Eq. (JTJ), becomes 

m = ^ < 1 . (2) 

^From this relation one obtains that for bombarding energies around and above 100 
MeV/nucleon, states with energy up to 10-20 MeV can be readily excited. The experi- 
mental problem is to ensure that the collision impact parameter is such that the nuclei do 
not overlap their matter distributions so that the process consists of Coulomb excitation 
only. This has been achieved by a careful filtering of the experimental events in terms 
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of scattering angles,. multiplicity of particles, angular distributions, using light and heavy 
targets, etc Q, 3, fl fl S H3- 

The theory of Coulomb excitation in low-energy collisions is very well understood |l| . It 
has been used and improved for over thirty years to infer electromagnetic properties of nuclei 
and has also been tested in the experiments to a high level of precision. A large number of 
small corrections are now well known in the theory and are necessary in order to analyze 
experiments on multiple excitation and reorientation effects. 

The standard semiclassical theory of Coulomb excitation at low energies assumes that 
the relative motion takes place on a classical Rutherford trajectory, as long as the transition 
energy Efi = Ef — Ei is small compared to the kinetic energy of the system. The cross 
section for exciting a definite final state | /) from the initial state | i) is then given by 



where Pi^f is the probability, evaluated in perturbation theory, o: 



the excitation of the 



target by the time-dependent electromagnetic field of the projectile 

In the case of relativistic heavy ion collisions pure Coulomb excitation may be distin- 
guished from the nuclear reactions by demanding extreme forward scattering or avoiding 
the collisions in which violent reactions take place [2j. The Coulomb excitation of relativis- 
tic heavy ions is thus characterized by straight-line trajectories with impact parameter b 
larger than the sum of the radii of the two colliding nuclei. A detailed calculation of rel- 
ativistic electromagnetic excitation on this basis was performed by Winther and Alder 
As in the non-relativistic case, they showed how one can separate the contributions of the 
several electric (EX) and magnetic (MA) multipolarities to the excitation. Later, it was 
shown that a quantum theory for relativistic Coulomb excitation leads to minor modifica- 
tions of the semiclassical results ^(|. In Ref. the connection between the semiclassical 

and quantum results was fully clarified. More recently, a coupled-channels description of 
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relativistic Coulomb excitation was developed |l2j]. 

The semiclassical theory of Coulomb excitation for low energy collisions accounts for 
the Rutherford bending of the trajectory, but relativistic retardation effects are neglected, 
while in the theory of relativistic Coulomb excitation recoil effects on the trajectory are 
neglected (one assumes straight-line motion), but retardation is handled correctly. In fact, 
the onset of retardation brings new important effects like the steady increase of the excitation 



cross sections with increasing bombarding energy. In a heavy ion collision around 100 
MeV/nucleon the Lorentz factor 7 is about 1.1. Since this factor enters in the excitation 
cross sections in many ways, as in the adiabacity parameter, Eq. (J2J), one expects that some 
sizeable (10 — 20 %) modification of the theory of nonrelativistic Coulomb excitation would 
occur. Also, recoil corrections are not negligible, and the relativistic calculations based on a 
straight-line parameterization of the trajectory are not completely appropriate to describe 
the excitation probabilities and cross sections. 

These questions are very relevant, as Coulomb excitation has proven to be a very useful 
tool in the investigation of rare isotopes in radioactive beam facilities [t|. Thus, it is appro- 
priate to investigate the effects of retardation and recoil corrections in Coulomb excitation 
at intermediate and high energies. In this article we will assess these problems by using the 
semiclassical approach of Ref . [13J . As we shall show in the next sections, both retardation 
and recoil effects must be included for bombarding energies in the range 30-200 MeV per 
nucleon. 

This can be accomplished in a straight-forward way in the semiclassical approach with 
a relativistic trajectory, appropriate for heavy ion collisions, and the full expansion of the 
electromagnetic propagator |l3j]. In most situations, the Coulomb excitation is a one-step 
process, which can be well described in first-order perturbation theory. Exceptions occur 
for very loosely bound nuclei, as for example the excitation of 11 Li Q, or 8 B 0, 3, m 
which case the electromagnetic transition matrix elements are very large due to the small 
binding and consequent large overlap with the continuum wavefunctions. Another exception 
is the excitation of multiple giant resonances, due to the strong collective response of heavy 
nuclei to the short electromagnetic pulse delivered in heavy ion collisions at relativistic 
energies 0, 0| . 

II. COULOMB EXCITATION FROM LOW TO HIGH ENERGIES 

In the semiclassical theory of Coulomb excitation the nuclei are assumed to follow classical 
trajectories and the excitation probabilities are calculated in time-dependent perturbation 
theory. At low energies one uses Rutherford trajectories while at relativistic energies one 

1m 



uses straight-lines for the relative motion 



In intermediate energy collisions, where 



one wants to account for recoil and retardation simultaneously, one should solve the general 
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classical problem of the motion of two relativistic charged particles. A detailed study of 
these effects has been done in Refs. Q, Q|. In Ref. Q] it was shown that the Rutherford 
trajectory is modified as the retardation due to the relativistic effects starts to set in. This 
already occurs for energies as low as 10 MeV/nucleon. It was also shown that if we use 
the scattering plane perpendicular to the z-axis, we may write the new Coulomb trajectory 
parameterized by 



x = a [coshu; + e] ; y = a Ve 2 — 1 sinhw; z = 0; t = — [w + e sinhu;] , (4) 



where e = 1/ sin (0/2), with O being the deflection angle. The impact parameter is related 
to the deflection angle by b = a cot (0/2). The only difference from Eq. (J1J and the usual pa- 
rameterization of the Rutherford trajectory at non-relativistic energies is the replacement of 
the half- distance of closest approach in a head-on collision ao = ZpZxe 2 jm^v 2 by a = ao/j. 
This simple modification agrees very well with numerical calculations based on the Darwin 
Lagrangian and the next order correction to the relativistic interaction of two charges Q|. 

Retardation also affects the dynamics of the Coulomb excitation mechanism and needs 
to be included in collisions with ener gies around 100 MeV/nucleon and higher. A detailed 
account of this has been given in Ref. [l3|. The end result is that the amplitude for Coulomb 
excitation of a target from the initial state | i) to the final state | /) by a projectile with 
charge Zp moving along a modified Rutherford trajectory is given by 

a /* = 2ATT (_1)M i S ( EX ^) M M EX i -ti) + S(M\,fi)M fi (M\, -//)} , (5) 

where A4fi(jrX, //) are the matrix elements for electromagnetic transitions, defined as 

M fi (EX,fi) = f^f^ J j/i(r) ■ V x L{ja(«t) Y^frfi} d 3 r (6) 
where L = — ir x V and 

Mfi(M\,fi) = J j/?(r) ' L ^ Kr ) ^,0)} d 3 r , (7) 

with lu defined as the excitation frequency Hu = Ef — Ei and k = lj/c. Using the Wigner- 
Eckart theorem 

M fl (7r\,-^ = (-l) 1 ^ I ' ; " I (//• ||..M(7rA)|| ./,) . (Si 




the geometric coefficients can be factorized in Eq. (J3J). 

The orbital integrals S(ir\, /x) are given by (after performing a translation of the integrand 
by w — > w + i (vr/2)) 



S(EX,n) 
S{M\n) 



% J(£A,/x), 



(9) 



Aca A 



[(2A + 1)/(2A + 3)] 1/2 [(A + l) 2 - ^ 2 ] 1/2 cot (#/2) /(MA, /x) 



with 



2A + 1 y/(A-/x)!(A + /x)! a+ttV2 
4tt (A-/x)!!(A + /i)!! 1 ; 




and 1 



I(E\,fJL)=-i(-± 



-tC/2 



c A(2A-1)!! 



for A + /i = even 
for A + /i = odd 



(e + z sinh w — -\/e 2 — 1 cosh 



w 



sinh w + 1)^ 1 



x 



(A + l)h\ — zh\ + \ e( coshw • h\ 

c 



and 



/(MA, ji) = ^ C/C) , A !! e"-'" 2 / <hc -'"V-"' 



(2A-1)!! 
(e + i sinh u> — -\/e 2 — 1 cosh 



w 



h x (z) 



(ie sinh w + 1)^ 

In the above equations, all the first-order Hankel functions h\ are functions of 



z = —C (ie sinh w + 1) 

c 



with 



lucl iua 

c = — 



(10) 



(ii) 



(12) 



(13) 



(14) 



v 7f 

The square modulus of Eq. (jSJ) gives the probability of exciting the target from the initial 
state | JjMi) to the final state | I/M/) in a collision with the center of mass scattering angle 



1 There is a misprint in the power of v/c in Eqs. 3.11 and 3.15 of Ref. [l^. 
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If the orientation of the initial state is not specified, the cross section for exciting the 
nuclear state of spin If is 

aV 1 



da, 



<->/ 



4 2L 



- I a /» I 2 dfi, (15) 

Mi,Mf 

where a 2 e 4 dQ/A is the elastic (Rutherford) cross section. Using the Wigner-Eckart theorem, 
Eq. (jHJ), and the orthogonality properties of the Clebsch-Gordan coefficients, gives 

^T = ^^ Ge ^ (2A + 1)3 |S(7rA ^ }l ' (16) 
where 7r = E or M stands for the electric or magnetic multipolarity, and 

B(-kX, h — ► If) — E I^M-^I 2 

= |(//||^(vrA)||/ 4 )| 2 (17) 

is the reduced transition probability. 



III. RELATIVISTIC AND NON-RELATIVISTIC LIMITS 

The non-relativistic limit is readily obtained by using v/c —>■ in the expressions in 
section II. In this case, z — > in Eq. (JT3J), and 

/ iA ^-i(2A-l)!!-l r 
(A + l)h x - zh x+1 - V -tC coshw • h x -> i\ (2A - , (18) 

which yields 

I(EX,fj,) = e"</ 2 r d„, e-Cecc-h^ (f + j smb™ = v^cosh^ 
J_ 00 (ze sinh w + 1) A +^ 

I(MX,fi) = I(EX + l,fi) . (20) 
These are indeed the orbital integrals for non-relativistic Coulomb excitation, as defined 



in Eq. (II-E.49) of Ref. |l5|. 



In the relativistic limit, v/c — > 1, £ — > in Eq. (|14j) and e ~ b/a — > oo. However, the 
combination 

£(&) = = ^ (21) 

7f 



can remain finite. 

The results for the orbital integrals can be expressed in closed analytical forms. First we 
translate back the integrands in Eqs. (fTT]) and (fT2j) by —i (ir/2) to get 

1 



1 !! 



(cosh w\ 



(A + l)h\ — zh\ + i + £—£ sinhw ■ h\ 

c 



and 



/(MA,//) 



i(v(/c 



,A+1 



(2A- 1)!! 



00 ^ fcA (z) e ie^ o+i^w; 

(cosher 



(22) 



(23) 



where now z — — £ coshw, and we took the limit £ — > and e — > 00. For the lowest multipo- 

c 

larities these integrals can be obtained in terms of modified Bessel functions by assuming the 
long- wavelength approximation, £ (77) <C 1, valid for almost all cases of practical interest. In 
this case, we can also use the approximation of Eqs. (|18|). From Eq. (j^J), in the relativistic 
limit, sinhw = vt/b and r = b cosh u>. Thus, the integrals can be rewritten as 



va x e-</ 2 



I(MX,fi) = va x+1 



dt e iivt/h 



{b + ivty 



(52 + w 2 t 2)(A+ M +l)/2 



.Vt V . 

Xb c 



/ 7te^ t/f> 



(b + k>i)^ 



'52 _|_ v 2 t 2)(A+ M +2)/2 • 

These integrals can be calculated analytically 1(| to give 



(24) 
(25) 



I(E\,fi) = F(A,/i,£) - —4 — . 



J(MA, yU ) = F(A + l,/i,0 (26) 



where 



2-1 * 



(&0 ^ 2" 



1 



p(n—n, n+fi) 
,ri—fi A— n 



(0) 



p(-M, A 6 ) (0) 



K n {t) 



(27) 



n=— A A 

In this equation P„ are the Jacobi polynomials, and K n (x) are modified Bessel functions. 
Since A + \i =even (odd) for electric (magnetic) excitations, we only need to calculate the 
integrals for /i = ±1 for the El multipolarity, /i = 0, ±2 for the E2 multipolarity, and \x = 
for the Ml multipolarity, respectively. 
To leading order in £, 



I(E1,±1) 



2a 



T(£), 7(M1,0) = 7(£2,0) 



2a' 



1(0, 7(E2,±2) = |!j(0 



x(0 = ^i(0 



for 
for 



£ < 1 
£ > 1 



(28) 
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When inserted in Eqs. (EJ) and (JTfij) the above results yield the correct relativistic Coulomb 
excitation cross sections [9j in the long wavelength approximation. Thus, we have shown ex- 
plicitly that the Equations (pf|)- (jl6)) reproduce the non- relativistic and relativistic Coulomb 
excitation expressions, as proved numerically in Ref. [13]. We can now analyze the interme- 
diate energy region (E^ab ~ 100 MeV/nucleon), where most experiments with radioactive 
beams are being performed. 



IV. GAMMA-RAY ANGULAR DISTRIBUTIONS 

As for the non-relativistic case 



171 ] , the angular distributions of gamma rays following 
the excitation depend on the frame of reference considered. It is often more convenient to 
express the angular distribution of the gamma rays in a coordinate system with the z- 
axis in the direction of the incident beam. This amounts in doing a transformation of the 
excitation amplitudes by means of the rotation functions D 3 m , m . The final result is identical 
to the Eqs. (II. A. 66- 77) of Ref. Q], with the non-relativistic orbital integrals replaced by 
Eqs. (fTT|) and (|12jh respectively. The angular distribution of the gamma rays emitted into 
solid angle f2 7 , as a function of the scattering angle of the projectile (0, $), is given by 

W (Q 7 ) = £ a x kK (6, $, C) A^Y kK (Q,) . (29) 

kK 

In our notation, the z-axis corresponds to the beam axis, and the a^ K (0, $, £) are given by 

a x kK (9, d>, C) = b x kK /b x m , (30) 
where, for electric excitations jl^ ]. 



2 f X X k\ ^ ( X X k 

C A (e,$,0 = - 7 =— £ (-if , , 

V2TT1 I i -i o / ±f K , U V «' 



x Fa, (|, 0) F v (|, 0) 7(EA, a*)/(SA, //)£>^ (f + ~, ~, *) • (31) 
In Eq. ()29)1 the coefficients A k are given by 

4 A) = F k (A, /<, //) F k (!, I', I 9 : If) &i A /' , (32) 
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where |A;| 2 is the intensity (in sec x ) of the 2'-pole radiation in the 7-transition from the 
excited state g to the state /. Explicitly, the /-pole conversion coefficient A; is given by 



A 



^H^'r^^^^ii-^ii^. (33) 



j [(2/ + ly.iyft \cJ 

with s(l) = I for electric (it = E) and s(l) = 1 + 1 for magnetic (tt = M) transitions. The 
product A/A;/ is always real since (— l) 11 ^ = (the parity). The coefficients F k (I, I', I g , If) are 
geometrical factors defined by 



F k (I, I', I g , I f ) = (-1) x l(2l + 1) (21' + 1) (21 f + 1) (2k + 1) 

(34) 



( I V k^ 


H 


I I' k } 






, J / J / j 



and 

F k (l,I t ,If) = F k (l,l,I l ,I f ). (35) 

The normalization of the coefficients a kK (Q,Q,() is such that <2q (©,$,£) = 1. Only 
terms with even k occur in Eq. (|29j). 

The total angular distribution of the gamma rays, which integrates over all scattering 
angles of the projectile, is given by 

W (ft,) = Yl a k (0 A[ X) P k (cos ft,) , (36) 
k 

where the z-axis corresponds to the beam axis and the statistical tensors are given by 

4 (0 = *>M , (37) 

where (for electric excitations) fl^ . 

Pi 2 ^A A &;\ ^ „ f X X k 

x F v (|, 0) >V , 0) rfe e J(£A, fi)I(EX, p!)Y kK (|, - + sin" 1 (1/e) 

(38) 

eo is the minimum value of the eccentricity, associated with the maximum scattering angle 
@o by eo = 1/ sin (O /2). One can show that the coefficients b k x are real, even if the orbital 
integrals are not. Their imaginary parts cancel out in the sum over fi/j,'. 
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For Ml excitations 



bf 11 




V2k + 1 

de e (e 2 - l) /(Ml, 0)/(Ml, 0)y fc0 ( ^ , 3 + sin " 1 (V<0) ■ (39) 

The normalization of the coefficients a£ (Q is such that (£) = 1. Again, only terms 
with even a£ (£) occur in Eq. (jHEJ)- 

In the case of Ml excitations Eq. (|39|) contains only k = fi = y! = and one gets 
a ¥ 1 (0 — 1) independent of (. Since for small ( the magnitude of a£ (() decreases ap- 
preciably with k we will only consider gamma ray emission after excitation through electric 
multipoles, in particular, the dependence of af 1 and af 2 on (. This dependence is very weak 
at energies E iab > 100 MeV/nucleon. In that case, one can use the approximate relations, 
Eq. (|28|) . for excitation energies such that ^ C 1. This condition is met for reactions with 
neutron-rich or proton-rich nuclei where the excitation energies involved are of the order of 
E x ~ 1 MeV. It is then straightforward to show that 

af 1 = 1, af 2 = -2, and af 2 = -0.25. (40) 

We thus come to the important conclusion that in high energy collisions and low excitation 
energies, E x ~ 1 MeV, the angular distribution of gamma-rays from decays after Coulomb 
excitation does not depend on the parameters ( or £. 

Although the cross sections for Ml and E2 excitations do not contain interference terms, 
the 7-decay of the excited state can contain an interference term with mixed E2 + Ml 
multipolaritites. The angular distribution of the gamma rays from the deexcitation of these 
states is given by [l^ 

W E2 ' M1 (0 7 ) = 2 V /^ V /^]T af > M1 (C) F k (1, 2, 4 I f ) ^Ai'Fk {I, I', I a , If) Pk (cos ft, 



'77 ) 

k IV 

(41) 

where gm\ {^et) is the total magnetic dipole (electric quadrupole) excitation cross section 
and where the sign of the square root is the same as the sign of the reduced matrix ele- 
ment (i (Ml)|| /) ((i (•£'2)11 /)). These latter are the same as those occurring in the 
radiative decay / — > i. The af 2 ' M1 coefficients in Eq. (J41J) are given by 



E2,M1 _ iE2,Ml 



vV^V^ ' (42) 
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where is given by Eq. (J35j) and 

^E2,M1 4 / 2 1 k \ 

fc V/2T+1 1-10 ^ 



i' 

oo 




X 



Y 2 » (~ 0) F 20 (|, 0) f° rfe e v 7 ^! I(E2, fi)I(Ml, 0)Y kfl (~ | + sin" 1 

^ CO 



(43) 

Using the relations in Eq. (}28|) and performing the summation above it is straightforward 
to show that for £(-R) -C 1 

E2,M1 n r i 

a k = (J, tor any k . 

Thus, in high-energy collisions, there is no interference term from mixed E2-M1 excitations 
in the angular distribution of emitted gamma rays. 

The form of the expressions for the angular distribution used here has followed that 
of Chapter 11 of Ref. [15]. Experimenters usually find it more convenient to write the 
angular distribution in a slightly different form which separates the statistical tensors that 
describe the orientation of the state due to the excitation process from the geometrical 
factors associated with the 7-ray decay and gives the geometrical factors the same form as 
occurs in the formulation of 7-7 correlations (cf. Ref. (2 P a g e 311; see also Refs. 18, 1^). 
The general expression for the 7-ray decay into solid angle f2 7 after projectile scattering to 
the angle (0, $), where the transition takes place between the Coulomb-excited state / and 
a lower state g (see Eq. (|2~H|) ) becomes: 

W (fi 7 ) = Y, a L (©, A k {6 7 ll'I g I f )Q k {E 7 )Y kK (fi 7 ) , (44) 

where the Ak^U'Iglj) coefficients are related to the so-called F-coefficient (Eq. (JSJ)) for 
the 7-ray transition between the states If and I g with mixed multipolarities / and /' and 



mixing ratio (L 
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by the expression 



A k (S 7 ll'I f I g ) = [F k (llI f I g ) + 25 7 F k (ll'I f I g ) + %F k (l'l'I f I g )]/(l + (45) 

Note that for k = we have A Q = F = 1, and due to the normalization used, the matrix 
elements, Eq. (J3l?j) . are not needed. In most situations one is interested in the possible 
mixing of E2 and Ml multipolarities in the decay of / — > g. Thus, one only needs the 
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E2/M1 mixing ratio. The quantity Qk(Ej) is the solid-angle attenuation coefficient which 



takes account of the finite solid-angle opening of the 7-ray detector 



2l|. It follows that 



a x kK (0, $, C) = a x kK (6, $, C) F k (\\lilf). (46) 

In the case where the particle scatters into an annular counter about the beam direction 
(i.e. the z-axis), or for angular distributions where the scattered particle is not detected at 
all (Eq. (jSEJ)), only k = terms survive. Usually the coefficients are normalized so that 

W(9 7 ) = l+J2 B k (0 Ak&ll'IglAQkiEjPk (cos# 7 ) . (47) 

fc=2,4 

The alignment of the initial state is now specified by the statistical tensor B^ which is related 
to the statistical tensors introduced above by 



Bl = V2k + 1 a x k0 /a x 00 , (48) 

when the particle is detected in an annular counter, and 

B x k =aUC)F k (\\IJ f ), (49) 

when the particle is not detected at all (or detected in such a way as to include all kinemat- 
ically allowed scattering angles). 

V. NUMERICAL RESULTS 

In Table H] we show the numerical results for the orbital integral I(E2,fi) for a deflec- 
tion angle of 10° and for ^ — 2,0, —2. The calculations have been done using the code 
. The results for 7 = 1 agree within 1/1000 with the numerical values ob- 

1 n 

23J, also reprinted in Table 11.12 of Ref. One observes that the results of 
the integrals for 7 = 1.1, corresponding to a laboratory energy of about 100 MeV/nucleon, 
differ substantially from the results for 7 = 1 (non-relativistic), specially for large values of 
(. For a fixed scattering angle ( increases with the excitation energy. Thus, one expects 
that the relativistic corrections are greater as the excitation energy increases. 

For 7 = 1 the imaginary part of the orbital integrals vanishes. But as ( and 7 increase the 
imaginary part becomes important. This is shown in Fig. Q]where the ratio of the imaginary 



COULINT |22 
tained in Ref. 
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to real parts of the orbital integral I(E2, 2) is shown for ( = 0.1 (dashed curve) and ( = 1 
(solid curve) as a function of 7. 

Except for the very low energies such that ao attains a large value, and for the very large 
excitation energies fiw, the parameter ( is much smaller than unity. Also, at intermediate 
energies the scattering angle is limited to very forward scattering. It is useful to compare the 
orbital integrals with their limiting expressions given by Eq. (|2*Hj) . i.e. the relativistic limit 
to leading order in £. This is shown in Fig. Ela) where the real (solid lines) and imaginary 
parts (long-dashed lines) of the orbital integral I (El, 1) are compared to the approximation 
of Eq. (dashed line) for 7 = 1.1 (Ei ah ~ 100 MeV/nucleon). Figure Efb) shows the 
same results, but for the orbital integral I(E2,2). The comparison is made in terms of the 
variable £ = eQ which is the appropriate variable for high energy collisions. Only for £ <C 1 
do the expressions in the relativistic limit reproduce the correct behavior of the orbital 
integrals. Although the imaginary parts of the orbital integrals are small, the real parts 
show substantial deviations from the approximations of Eq. (128)1 at intermediate energies 
(E lah ~ 100 MeV/nucleon). 

We now apply the formalism to specific cases. We study the effects of relativistic cor- 
rections in the collision of the radioactive nuclei 38 ' 40 > 42 S and 44 ' 46 Ar on gold targets. These 
reactions have been studied at Ei a b ~ 40 MeV/nucleon at the MSU facility [24] . In the fol- 
lowing calculations the conditions may be such that there will be contributions from nuclear 
excitation, but these will be neglected as we only are interested in the relativistic effects 
in Coulomb excitation at intermediate energy collisions. In Table ITT1 we show the Coulomb 
excitation cross sections of the first excited state in each nucleus as a function of the bom- 
barding energy per nucleon. The cross sections are given in milibarns. The numbers inside 
parenthesis and brackets were obtained with pure non-relativistic and relativistic calcula- 
tions respectively. The minimum impact parameter is chosen so that the distance of closest 
approach corresponds to the sum of the nuclear radii in a collision following a Rutherford 
trajectory. One observes that at 10 MeV/nucleon the relativistic corrections are important 
only at the level of 1%. At 500 MeV/nucleon, the correct treatment of the recoil corrections 
(included in the equations (jllj) and ([12)1 ) is relevant on the level of 1%. Thus the non- 
relativistic treatment of Coulomb excitation Il7|] can be safely used for energies below about 

fl 

10 MeV/nucleon and the relativistic treatment with a straight-line trajectory [9j is adequate 
above about 500 MeV/nucleon. However at energies around 50 to 100 MeV/nucleon, accel- 
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c 


\,H = 2,2 


\,H = 2,0 


X,n = 2,-2 


0.0 


5.064(-3) [5.052(-3)] 


1.332(-2) [1.332(-2)] 


5.064(-3) [5.073(-3)] 


0.1 


8.675(-4) [1.105(-3)] 


6.505(-3) [7.621(-3)] 


1.195(-2) [1.205(-2)] 


0.2 


1.895(-4) [2.897(-4)] 


2.280(-3) [3.276(-3)] 


7.765(-3) [9.487(-3)] 


0.3 


4.425(-5) [7.849(-5)] 


7.311(-4) [1.296(-3)] 


2.245(-3) [5.470(-3)] 


0.4 


1.069(-5) [2.122(-5)] 


2.245(-4) [4.920(-4)] 


1.468(-3) [2.728(-3)] 


0.5 


2.637(-6) [5.599(-6)] 


6.716(-5) [1.821(-4)] 


5.442(-4) [1.251(-3)] 


0.6 


6.598(-7) [1.402(-6)] 


1.975(-5) [6.627(-5)] 


1.908(-4) [5.431(-4)] 


0.7 


1.668(-7) [3.169(-7)] 


5.739(-6) [2.383(-5)] 


6.438(-5) [2.268(-4)] 


0.8 


4.250(-8) [5.580(-8)] 


1.652(-6) [8.492(-6)] 


2.110(-5) [9.201(-5)] 


0.9 


1.090(-8) [1.906(-9)] 


4.724(-7) [3.004(-6)] 


6.765(-6) [3.650(-5)] 


1.0 


2.807(-9) [-5.003(-9)] 


1.343(-7) [1.057(-6)] 


2.131(-6) [1.422(-5)] 


1.2 


1.886(-10) [-1.845(-9)] 


1.071(-8) [1.291(-7)] 


2.031(-7) [2.074(-6)] 


1.4 


1.282(-11) [-3.551(-10)] 


8.426(-10) [1.556(-8)] 


1.858(-8) [2.902(-7)] 


1.6 


8.788(-13) [-5.699(-ll)] 


6.566(-ll) [1.857(-9)] 


1.651(-9) [3.941(-8)] 


1.8 


6.068(-14) [-8.381(-12)] 


5.078(-12) [2.200(-10)] 


1.433(-10) [5.227(-9)] 


2.0 


4.213(-15) [-1.170(-12)] 


3.904(-13) [2.591(-11)] 


1.222(-11) [6.809(-10)] 


4.0 


1.294(-26) [-1.211(-21)] 


2.362(-24) [l.lll(-20)] 


1.464(-22) [5.663(-19)] 



TABLE I: The classical orbital integrals for E2 Coulomb excitation. The Table lists the values 
of the classical orbital integrals I(E2, fx) for a deflection angle of 10° and for ^ = 2,0, —2. These 
entries are given in the form of a number followed by the power of ten which it should be multiplied. 
The value outside (inside) the brackets are for 7 = 1 (7 = 1.1). 

erator energies common to most radioactive beam facilities (MSU, RIKEN, GSI, GANIL), it 
is very important to use a correct treatment of recoil and relativistic effects, both kinemat- 
icaliy and dynamically. At these energies, the corrections can add up to 50%. These effects 
were also shown in Ref. for the case of excitation of giant resonances in collisions at 
intermediate energies. As shown here, they are also relevant for the low-lying excited states. 

As another example, we calculate the Coulomb excitation cross sections of n Be projectiles 
on lead targets. 11 Be is a one neutron halo- nucleus with one excited bound state (| at 
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CD 




FIG. 1: The ratio of the imaginary to real parts of the orbital integral I(E2, 2) is shown for ? = 0.1 
(dashed curve) and ( = 1 (solid curve) as a function of 7. 

320 keV). Its | + ground state is strongly coupled to the excited state with the strongest 
El transition observed between bound nuclear states. The B(E1) value for this transition 
is 0.116 e 2 fm 2 25]. The one-neutron separation energy is only 506 keV and the coupling 
to the continuum has to be included in an accurate calculation. However the influence of 
higher-order effects in the Coulomb excitation of the excited state in intermediate energy 
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271 128| to be less than 7%. We thus neglect these effects 



collisions was shown in Refs. 
here. 

In Ref. a recoil correction for the theory of relativistic Coulomb excitation was pro- 
posed. It was shown that one can use the equations for relativistic Coulomb excitation and 
obtain reasonable results for collisions at low energies if one replaces the impact parameter 
b by 

b' = b+^-a. (50) 

The advantage of this approximation is that one can use the analytical formulas for relativis- 
tic Coulomb excitation (e.g., Eqs. EBKZt and easily include the recoil correction Eq. (J5U|) . 
We define the percent deviation 

('-'"exact ®i 



A,; 



C^cxact 



(51) 
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where <j{ is the cross section obtained with the relativistic (i = R), non-relativistic (i = NR), 
and (i = RR) with the relativistic equations for Coulomb excitation Eq. (ffilj) but with the 
recoil correction Eq. (J5U|) . respectively. Figure 121 is a plot of Eq. (|51|) for the excitation 
of the 0.89 MeV state in 40 S + 197 Au collisions clS db function of the bombarding energy. 
One observes that the largest discrepancy is obtained by using the non-relativistic equations 
(NR) for the Coulomb excitation cross sections at high energies. The relativistic analytical 
equations (R) also do not do a good job at low energies, as expected. But the relativistic 
equations with the recoil correction of Eq. ()50|) improve considerably the agreement with the 
exact calculation. At 10 MeV/nucleon the deviation from the exact calculation amounts to 
6% for the case shown in Fig.|21 However, the deviation of the RR treatment tends to increase 
for cases where higher nuclear excitation energies are involved ^J. The cross section for 
the excitation of low energy states is mainly due to collisions with large impact parameters 
for which recoil corrections are not relevant. For high lying states, e.g. giant resonances, 
only the smaller impact parameters are effective in the excitation process. Therefore, in this 
situation, the correct treatment of recoil effects is more relevant. 

For 11 Be projectiles on lead targets at 50 MeV/nucleon the Coulomb excitation cross 
sections of the excited state in n Be are given by 311 mb, 305 mb and 398 mb for non- 
relativistic, exact, and relativistic calculations, respectively. At 100 MeV/nucleon the same 
calculations lead to 159 mb, 185 mb and 225 mb, respectively. Thus, the same trend as in 
the results of Table |H] is also observed for El excitations of low-lying states. 

Experiments of Coulomb excitation of n Be projectiles have been performed at 
GANIL (43 MeV/nucleon) Q, at RIKEN (64 MeV/nucleon) H and at MSU (57-60 
MeV/nucleon) 31]. The extracted values of B(E1) ~ 0.05 e 2 fm 2 in the GANIL experi- 
ment is in disagreement with the lifetime experiment of Ref. 2^ and could not be explained 
by higher-order effects in Coulomb excitation at intermediate energy collisions j3, 27. I^. 



However, the deduced values of B(E1) ~ 0.1 e 2 fm 2 in the RIKEN and MSU experiments 

25l and with the theoretical cross 



are in good agreement with the lifetime measurement 
sections of Coulomb excitation. 

We now study the effects of retardation in the angular distributions of gamma-ray decay- 
ing from Coulomb excited states. We first test the range of validity of the approximation in 
Eq. ([40)1 . For this purpose we artificially vary the energy of the first excited 2 + state in 38 S. 
This would simulate what happens in the case of a nucleus with very low-lying excited states, 
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E lab [MeV/nucleon] 



FIG. 2: Eq. (j5D for the excitation of the 0.89 MeV state in 40 S + 197 Au collisions as a function 
of the bombarding energy. The solid line corresponds to the use of the non-relativistic integrals 
Eq. (|19|) compared to the exact calculation using Eq. The same is plotted for the other two 

cases: (R) with the relativistic Eq. (|26|) . and (RR) with the relativistic Eq. (|26|) using the recoil 
correction Eq. 1)50(1 . 

or very high bombarding energies. As we see from figure 0J the statistical tensors, B^ 2 ^ , 
asymptotically attain the values B^ 2 ^ = 1.19 and £?j = 0.267 according to the limits in 
Eq. flUJ) and the definition of B\ in Eq. flU}, since F 2 = -0.5976 and F 4 = -1.0690. The 
convergence to these asymptotic values increases with the bombarding energies, as expected 
from the conditions which lead to validity of EJj/^hc <C 1, for the lowest impact parameters 
b, which are the most relevant for the Coulomb excitation process. At 1 GeV/nucleon this 
condition is easily met for states of the order of 1 MeV. 

Finally, we show in Table lTTTl the statistical tensors Bf 2 in equation (jTTj) for 38 S projectiles 
at 100 MeV/nucleon incident on gold targets and scattering to all kinematically allowed 
angles with closest approach distance larger than the sum of the nuclear radii. NR (R) 
denotes the non-relativistic (relativistic) values. We notice that the statistical tensors are 
not as much influenced by the retardation and recoil corrections as in the case of the cross 
sections. The reason is that the statistical tensors involve ratios of the integral of the orbital 
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FIG. 3: Upper figure: Real (solid lines) and imaginary part (long-dashed lines) of the orbital 
integral I (El, 1) for 7 = 1.1 (-Ei a b — 100 MeV/nucleon). The approximation of Eq. (|28[) is shown 
by the dashed line. Lower figure: Same plot, but for the orbital integral I(E2,2). 

integrals. These ratios tend to wash out the corrections in the orbital integrals due to 
relativistic effects. On the other hand, for scattering to a specific angle the corrections can 
be larger because the a^ K in Eq. (jUj) are determined largely by geometry and hence they 
can be sensitive to both relativistic distortions of the orbit and recoil effects causing non 
straight-line trajectories. 



VI. CONCLUSIONS 

We have extended the study of Ref. [k| to include retardation effects in the Coulomb 
excitation of low-lying states in collisions of rare isotopes at intermediate energies (E^ ~ 100 
MeV/nucleon). In particular, we have studied the effects of retardation and recoil in the 
orbital integrals entering the calculation of Coulomb excitation amplitudes. We have shown 
that the non-relativistic and relativistic theories of Coulomb excitation are reproduced in 
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Nucleus 


Ex 

[MeV] 


B(E2) 
[e 2 fm 4 ] 


10 MeV/ A 
ac [mb] 


50 MeV/ A 
ac [mb] 


100 MeV/ A 
ac [mb] 


500 MeV/ A 
ac [mb] 


38g 


1.29 


235 


(492) 500 [651] 


(80.9) 91.7 [117] 


(40.5) 50.1 [57.1] 


(9.8) 16.2 [16.3] 


40g 


0.89 


334 


(877) 883 [1015] 


(145.3) 162 [183] 


(76.1) 85.5 [93.4] 


(9.5) 20.9 [21.] 


42 S 


0.89 


397 


(903) 908 [1235] 


(142.7) 158 [175] 


(65.1) 80.1 [89.4] 


(9.9) 23.2 [23.4] 


44 Ar 


1.14 


345 


(747) 752 [985] 


(133) 141 [164] 


(63.3) 71.7 [80.5] 


(8.6) 17.5 [17.6] 


46 Ar 


1.55 


196 


(404) 408 [521] 


(65.8) 74.4 [88.5] 


(30.2) 37.4 [41.7] 


(5.72) 10.8 [11] 



TABLE II: Coulomb excitation cross sections of the first excited state in 38 > 40 > 42 S and 44 > 46 Ar projec- 
tiles at 10, 50 100 and 500 MeV/nucleon incident on gold targets. The numbers inside parenthesis 
and brackets were obtained with pure non-relativistic and straight-line relativistic calculations, 
respectively. The numbers at the center are obtained with the full integration of equations (jll|) 
and (H2J). 





NR 


Exact 


R 


BP 


0.95 


1.03 


1.11 


BP 


0.183 


0.192 


0.207 



TABLE III: Statistical coefficients entering equation (jlTj) for 38 S projectiles at 100 MeV/nucleon 
incident on gold targets. 

the appropriate energy regime. We have also shown that at intermediate energies corrections 
to the low- or high-energy theories of Coulomb excitation are as large as 20%. 

We have studied the excitation of the first excited states in n Be, 38 ' 40 ' 42 S and 44,46 Ar 
projectiles incident on gold and lead targets. It is clear from the results that retardation 
corrections are of the order of 10%-20% at bombarding energies around 50-100 MeV/nucleon. 
Therefore, they must be accounted for in order to correctly analyze the cross sections and 
angular distributions of decaying gamma-rays in experiments at radioactive beam facilities 
running at intermediate energies. 

Another important consequence of our study is that retardation effects must also be 
included in calculations of higher-order effects (e.g., coupled-channels calculations), common 
in the Coulomb breakup of halo nuclei [4]. Work in this direction is in progress. 
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FIG. 4: Statistical tensors Bf> and B\ ' as a function of the excitation energy of the lowest 
state in 38 S projectiles incident on gold targets. The energy of the state is varied artificially. 
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